Effect of interactions on harmonically confined Bose- Fermi mixtures in optical lattices 
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We investigate a Bose-Fermi mixture in a three-dimensional optical lattice, trapped in a har- 
monic potential. Using Generalized Dynamical Mean-Field theory, which treats the Bose-Bose and 
Bose-Fermi interaction in a fully non-perturbative way, we show that for experimentally relevant 
parameters a peak in the condensate fraction close to the point of vanishing Bose-Fermi interaction 
is reproduced within a single band framework. We identify two physical mechanisms contributing 
to this effect: the spatial redistribution of particles when the interspecies interaction is changed 
and the reduced phase space for strong interactions, which results in a higher temperature at fixed 
entropy. 
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Cold atomic gases offer the exciting possibility to re- 
alize mixtures of fermions and bosons in optical lattices 
[US]. This yields a very interesting system without a 
clear analog in conventional solid state physics. The 
depth of the optical lattice and the Bose-Fermi inter- 
action can be tuned independently, which has made it 
possible to explore this system experimentally in a de- 
tailed way [4 . One of the key questions that has been 
investigated is the effect of the fermions on the mobil- 
ity of the bosons. When the fermions are slow compared 
to the bosons, they act as (dynamical) impurities. Fast 
fermions, on the other hand, mediate a long-range inter- 
action between the bosons. In both cases the fermions 
induce a shift of the bosonic superfluid-Mott transition, 
which has been intensively studied by different experi- 
mental groups [1, 3, 4 . Most importantly, these time-of- 
flight experiments showed significant loss of bosonic co- 
herence when fermions are present, thus indicating that 
adding the fermions stabilizes the Mott insulator phase 
[1, 3 . A pronounced peak in the bosonic visibility as a 
function of the Bose-Fermi interaction was observed close 
to vanishing interaction [4] . Theoretical studies trying to 
explain this phenomenon [5HT3]. however, predict that 
within a single-band approximation at zero temperature 
the Mott region shrinks when fermions are added. Heat- 
ing is proposed as an explanation of this apparent contra- 
diction [13] and also multiband-effects can lead to an ex- 
tension of the Mott insulating regime (in a homogeneous 
system) , due to the self-trapping effect [TQtiT2| [T4] . How- 
ever, in the experiment reported in [4] the increase of vis- 
ibility happens for small Bose-Fermi interaction, where 
multi-band effects are supposed to be negligible. 

In this Letter we investigate this problem by means 
of Generalized Dynamical-Mean Field Theory (GDMFT) 
[T5] . This is a non-perturbative method which is well- 
controlled in three dimensions. We will show that in 
the homogeneous limit in the single-band framework we 
indeed reproduce the shrinking of the Mott lobes cor- 



responding to an induced attractive interaction between 
the bosons mediated by the fermions. To reproduce the 
experimental results we include two essential physical in- 
gredients. First of all, we take into account the (experi- 
mentally always present) confining potential, which leads 
to an inhomogeneous system. This allows us to probe 
the interaction-dependent distribution of particles in the 
trap and its effect on the bosonic coherence. Secondly, 
we account for the fact that the experiments are not per- 
formed at constant temperature, but at constant entropy 
(assuming that the ramp- up of the optical lattice is adi- 
abatic). For strong attractive/repulsive interactions the 
phase space is strongly reduced, which results in a lower 
entropy at constant temperature and hence a higher tem- 
perature at constant entropy. We will indeed show that 
the temperature has a minimum close to the position of 
the peak in the condensate fraction. A combination of 
these two effects is crucial for a qualitative understanding 
of the experimental observations. 

We study a Bose-Fermi mixture within the single-band 
approximation and the tight-binding limit, described by 
the Hamiltonian 

H = -E{*/ c i c i+^'} + ^ (V ai -na)nf 

(ij) i,a=bj 
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where c\ (b\) is the fermionic (bosonic) creation opera- 
tor at site z, while n{ = cJq (n\ = b\bi) denotes the 
corresponding number operator. Hf(b) is the chemical 
potential for fermions (bosons). C/j, and Ubf are the on- 
site boson-boson and boson-fermion interaction parame- 
ters, which are related to the bosonic scattering length a& 
and the Bose-Fermi scattering length a^f by the overlap 
integral of the corresponding Wannier functions. The 
summation index denotes summation over near- 

est neighbors, and tf(p\ is the tunneling amplitude for 
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fermions (bosons). 

The harmonic trapping potential is incorporated 
through the term V ai = V a \ri — r a o\ 2 (a = b,f). Here 
we assume the positions of the minima to be equal 
(rjo = 1*60 )? thus ignoring a possible differential sag be- 
tween the two clouds. The parameters we have chosen in 
the simulations for trapped mixtures are closely related 
to the experimental ones in Ref. [4 . The two species 
in this experiment are bosonic 87 Rb and fermionic 40 K. 
By working at an optical lattice wavelength A = 755 
nm, it is ascertained to have equal dimensionless lat- 
tice depth for both species V&jE b R = V^jE f R (Vg L 
being the optical lattice depth and E R the recoil en- 
ergy). This fixes the ratio between the hopping con- 
stants as tfjth — m^/ruf = 2.2. The intrabosonic repul- 
sion is fixed by the bosonic background scattering length 
of a\) ~ lOOao (ao being the Bohr radius), whereas the 
interspecies scattering length is tunable by a Feshbach 
resonance. We fix the total number of particles to ap- 
proximately Nb = 4 x 10 5 and iVj = 3x 10 5 , since these 
numbers yielded the most pronounced experimental sig- 
nal. The experimental trap is pancake-shaped. Since 
the dipole trap is counter-acted by the laser profile of 
the blue-detuned lattice beam, the precise trapping fre- 
quencies depend on the lattice depth. Our simulations 
are performed in an isotropic three-dimensional trap with 
trapping frequencies equal to the averaged experimental 
frequencies uj a = y {uj^_) 2 uj^. 

The single band approximation made here cannot be 
justified for the entire range of interactions realized in the 
experiments. Significant multi-band effects show up for 
large attraction between bosons and fermions [lOj [TT] . 
However, we focus on the regime of small Bose-Fermi 
interactions, where multi-band effects are supposed to 
be weak. In this way we are able to isolate the effects 
of the temperature and the spatial redistribution of the 
particles and to distinguish them from possible multi- 
band effects. 

To solve this interacting many-body system, we use 
GDMFT, which is a combination of Gutzwiller mean- 
field theory for the bosons and Dynamical Mean-Field 
Theory for the fermions [15 . The effect of the har- 
monic trap is treated within local density approximation 
(LDA): the harmonic potential is modelled by a spatially 
varying chemical potential. Exact diagonalization (ED) 
of the Anderson Hamiltonian is used as the DMFT im- 
purity solver. To calculate the entropy we employ the 
thermodynamic relation dVt = —SdT — Nfd/if — N^dji^ 
for the grand potential Q. 

Before turning to the trapped mixture, we first present 
GDMFT results for a homogeneous system: for constant 
fermionic filling rif = 0.7 the phase border between the 
superfluid and the Mott phase with bosonic filling = 
1 is shown in Fig. [I] for different values of the Bose- 
Fermi interaction. We observe that both for repulsive 
and attractive interaction the Mott lobe shrinks, thus 
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FIG. 1: (Color online) Phase border between bosonic Mott 
insulator and superfluid with rib = 1 for fermionic filling 
rif = 0.7, hopping ratio tf/tb = 2.2, temperature T — >• and 
various values oiUbf- The Hartree shift due to the Bose-Fermi 
interaction is subtracted from the bosonic chemical potential 
to facilitate an easy comparison. 

giving evidence that the GDMFT scheme fully includes 
the dynamical screening of the bosonic interaction by the 
fermions. 

We now turn to the inhomogeneous case, for which we 
focus on the condensate fraction of the bosons, which 
within LDA is obtained as 

N c = J d 3 r\(b(r)}\ 2 / J d 3 rn b (r). (2) 

The experimentally measured visibility is a (monotonic) 
function of N c within GDMFT, therefore a maximum in 
the visibility corresponds to a peak in N c . 

The condensate fraction is calculated at constant en- 
tropy per particle s = Stot/^V"tot- ^n estimate for the 
entropy per particle in the experiment is s = 0.86/cb, but 
since this value is not known very accurately, we consider 
a range of values for s. 

We focus on the case of Vq L /E r = 9 (shown in Fig. 
|2j, but data for VgjE% = 5 and Vg L /E% = 12 (see 
Fig. [3]) show similar features. Most importantly, the 
condensate fraction (Figs. [2^i, |3|i,b) shows a clear peak 
close to vanishing Bose-Fermi interaction for all lattice 
depths and for all values of the entropy considered, in 
agreement with the experimental observations [4 . 

To get insight into the physical mechanism behind the 
peak in 7V C , we plot the temperature at constant entropy 
in Fig. [2^ (cf. |3]3,d). Generally, the temperature has 
a minimum close to the position of the maximum in 
the condensate fraction Af c . In a qualitative way this 
minimum in the temperature at constant entropy can be 
understood as follows: for strong attraction only states 
in which bosons and fermions occupy the same lattice 
sites are energetically allowed. This means that the size 
of the accessible Hilbert space is strongly reduced and 
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FIG. 2: (Color online) a Total condensate fraction N c , e tem- 
perature T, and f spatial overlap of the bosonic and fermionic 
clouds as a function of Bose- Fermi scattering length a b f for 
various values of the entropy per particle and Vql/^S = 9 
(U b /tf = 5.2, Vf/tf = 0.0058, Vb/tf = 0.0064). b-d: Spa- 
tial profiles of the bosonic and fermonic density and the local 
condensate fraction N c (r) = \{b(r))\ 2 /Nb(r) as a function of 
distance from the trap center for three values of the Bose- 
Fermi interaction. 



hence the entropy at constant temperature is lowered. 
To obtain the same entropy, the temperature needs to 
be increased. The same mechanism works at strong re- 
pulsive interactions: in this case it is energetically unfa- 
vorable that a fermion and one or more bosons occupy 
the same site. This is an example of adiabatic heating if 
the Bose-Fermi interaction is changed. A similar effect 
has recently been observed in a fermionic mixture, where 
for increasing attractive interaction, the cloud expands 
due to a redistribution of entropy from spin to orbital 
degrees of freedom [16]. In the Bose-Fermi mixture, the 
effect is larger on the attractive than on the repulsive 
side, leading to a higher condensate fraction for repulsive 
than for attractive interactions, in agreement with the 
experimental findings [4]. This is because on the repul- 
sive side the bosons and fermions phase separate, whereas 
for attractive interactions all the particles are packed in 
the center. This spatial redistribution is shown in the 
radial density distributions in Fig. [2}>d where the den- 
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FIG. 3: (Color online) Condensate fraction (top row), temper- 
ature (middle row) , and overlap of the bosonic and fermionic 
clouds (bottom row) for Vq Li /E% = 5 (left column) and 
VqIj/Er = 12 (right column). Specific parameters are: 
Vol/^S = 5: Ub/tf = 1.06, Vf/tf = Vb/tf = 0.0024; 
VSJE% = 12: Ub/tf = 13.4, Vf/tf = 0.0108, V b /t f = 0.012. 



sity profiles for the bosonic and fermionic atoms and the 
local condensate fraction are plotted as function of the 
radius for strong attractive (b), weak repulsive (c), and 
strong repulsive (d) Bose-Fermi interaction. The profiles 
are plotted here at respective temperatures correspond- 
ing to total entropy equal to s = 0.86fc#. For strong 
Bose-Fermi attraction the bosons and fermions both oc- 
cupy the center of the trap. The fermions form a band 
insulator; the bosons have a high filling and are super- 
fluid. This implies that the bosons and fermions do not 
form bound states. If that were the case, the resulting 
composite particle would be fermionic in nature and the 
superfluid order parameter would vanish. For strong re- 
pulsion, on the other hand, the bosons occupy the center 
of the trap, whereas the fermions are pushed towards the 
outside. 

The most interesting regime is that of weak Bose-Fermi 
interaction, as shown in Fig. [2^. The value of a^f cho- 
sen here corresponds to the peak in N c . We observe that 
in this case the spatial distribution of the fermions is 
relatively fiat. This is due to the fact that for these pa- 
rameters the interparticle repulsion partly compensates 
the harmonic confinement, which in turn leads to a large 
phase space and the possibility to store entropy efficiently 
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as configurational entropy in the wings of the system. 
We can carry this analysis one step further by observing 
that this value of a& f is approximately distinguishing the 
regime of phase separation of the bosonic and fermionic 
clouds (for stronger repulsion) from overlapping bosonic 
and fermionic clouds (for weaker repulsion and attrac- 
tion). This is made clear in Fig. [2]p (cf. Fig. [3^,f), 
where the spatial overlap of the bosonic and fermionic 
clouds (defined by the integral J d 3 r Nb(r)Nf(r), here 
normalized to the total number of bosons Nt>) is plotted 
at constant entropy. The sharp decrease in overlap ap- 
proximately coincides with the minimum in temperature 
and the maximum in N c . It is indeed reasonable to argue 
that at the onset of phase separation the phase space is 
largest, since then the bosons and fermions are neither 
bound to the same lattice sites, nor to excluding sites. 
This leads to a low temperature at fixed entropy, which 
in turn causes a high condensate fraction. We find the 
maximum in N c , minimum in temperature, and onset 
of phase separation to happen at almost the same value 
of the interspecies interaction. However, the three do 
not exactly coincide, because of the complex interplay of 
many-body effects (leading also to a structure in 7V C at 
constant temperature when the interaction is changed) 
and the adiabatic heating effect explained above. 

In conclusion, we have investigated a trapped Bose- 
Fermi mixture by means of GDMFT. The dynamical 
screening of the Bose-Bose interaction by the fermions 
leads to decreasing Mott insulating lobes in the ground 
state of homogeneous mixtures. In the trapped case, a 
peak in the condensate fraction is obtained close to van- 
ishing Bose-Fermi interaction. We have shown that the 
dramatic decrease in condensate fraction for attractive 
Bose-Fermi interactions is caused by adiabatic heating 
due to the strong reduction of the available Hilbert space. 
For both positive and negative interspecies interaction 
strength we find a strong redistribution in the density 
of the two-component gas. The combination of both ef- 
fects leads to the peak in the condensate fraction close to 
vanishing Bose-Fermi interactions. We remark that these 
effects are universal and important for the interpretation 
of a large class of experiments with ultracold atoms in 
optical lattices. 

Although we thus obtained good qualitative agreement 
with the experimental results, there is no full quantita- 
tive agreement. In particular, it was observed that the 
peak in the visibility occurs exactly at vanishing Bose- 
Fermi interaction, whereas we find it to happen at a 
small, but positive value of the interaction. There are 
a number of possible explanations for this disagreement. 
First of all, our analysis relies on the assumption of adi- 
abatic loading of the optical lattice. Even for moderate 
interactions, however, one can expect mass transport to 
be affected, leading to longer adiabaticity timescales, as 
evidenced in recent transport experiments [T7J [18] . The 
strong particle redistribution for large attractive and re- 



pulsive interaction can easily induce additional heating 
and shift the peak in the condensate fraction. A second 
effect not included in our study is the presence of grav- 
itational sag between the bosonic and fermionic cloud 
in the experiment, due to their different masses. This 
separates the centers of mass of the two clouds already 
without interspecies interactions and additional repulsion 
enhances this separation. Since we observe in our calcu- 
lations that the point where the clouds separate is corre- 
lated with the maximum in the condensate fraction, the 
experimental sag is likely to shift the peak in the vis- 
ibility towards smaller Bose-Fermi interaction. Finally, 
residual uncertainties in the scattering length or multi- 
band effects could play a role in the exact location of the 
peak in the condensate fraction. Although large qual- 
itative effects are only expected for strong interspecies 
interactions [TQHT21 H4]. at small values of a^f multiband 
effects can still lead to renormalization of the parameters 
and small quantitative changes in the predictions. 
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Note added.- During completion of this work, related 
results were reported in [19]. 
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